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CN ■ Abstract 

I [ Magnetic Raman scattering from a frustrated spin-1/2 Heisenb erg-chain is 

U: 

t^ , considered with a focus on the uniform phase of the spin-Peierls compound 

CuGeOa- The Raman intensity is analyzed in terms of a Loudon-Fleury scat- 

'^ I tering process using a spinless-fermion mean-field theory developed for the 

O I frustrated spin-chain. A comparison to experimental data is presented and 
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the frustration and temperature dependence is studied. In good agreement 
with observed spectra a broad inelastic four-spinon continuum is found at 



fSJ ! low temperatures above the spin-Peierls transition. At high temperatures 

(N 

^^ I the intensity develops a quasi-elastic line analogous to experiment. 
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I. INTRODUCTION 

Elementary excitations of one-dimensional (ID) quantum spin-chains exhibit a num- 
ber of remarkable subtleties, such as the magnon continuum of the isotropic ID spin-1/2 
Heisenberg model or the excitation-gap for integer-spin systems as conjectured by Hal- 
dane. Recent and extensive studies of 'spin-ladder' compounds as well as of the germanate 
CuGeOs [|I|] have lead to renewed interest in quantum spin-chains. In this context, magnetic 
Raman scattering, if symmetry-allowed, has developed into an important tool to investigate 
the local spin dynamics [0-^. Here I will establish a simple framework to interpret the 
magnetic Raman scattering from a frustrated quantum spin-chain with a particular focus 
on the uniform phase of CuGeOs. 

CuGeOs is a quasi-lD anorganic spin-Peierls compound with a dimerization transition 
at a temperature Tsp — 14K [|l],0,||]. Its structure comprises of weakly coupled Cu02 chains, 
with copper in a spin-1/2 state 0. Since the nearest-neighbor (n.n.) exchange-coupling 
between copper spins along the Cu02 chains is strongly reduced by almost orthogonal in- 
termediate oxygen states |^ the next-nearest-neighbor (n.n.n.) exchange is relevant. Both, 
n.n. and n.n.n. exchange, are antiferromagnetic implying intra-chain frustration. A 
minimal model of CuGeOs is the J1-J2-S model 

H = J^J2[{l + {~y5)SrSl+, + aSl■Sl+2] • (1) 



Here, S; is a spin-1/2 operator, Ji ^ IQOK |T0|-[l2| is the n.n. exchange-coupling constant 



and 6 resembles the lattice dimerization which is finite for T < Tsp only, a is the intra- 
chain frustration-ratio a = J2/J1 where J2 is the n.n.n. exchange-coupling constant. A 
final consensus on the precise magnitude of a is still lacking. Studies of the magnetic 



susceptibility have resulted in a ^ 0.24 |jTO| as well as in a ?« 0.35 |TT[. The latter is 



consistent with a comparative investigation of magnetic susceptibility and thermal expansion 



12|. Therefore, very likely, CuGeOa displays a frustration induced contribution to the spin- 



gap at zero temperature irrespective of the actual lattice dimerization W^ . Depending on 



the magnitude of the frustration the zero-temperature dimerization S(T = 0) has to be on 
the order of 5(0) = 0.01. ..0.03 such as to enforce the size of the spin-gap observed in inelastic 
neutron scattering (INS) 0. 

Magnetic excitations in CuGeOs are clearly distinct among the uniform, i.e. T > Tsp, 
and the dimerized, i.e. T < Tsp, phase. While the dynamic structure factor exhibits a 



gapless two-spinon continuum similar to that of the ID Heisenberg chain above Tsp []14 



well defined magnon-like excitations have been observed below Tsp [|r^]. These magnons 
reside within the spin-gap and are split off from the two-spinon continuum. They have been 



interpreted as two-spinon bound states |[T6|JT7|1 



Magnetic Raman scattering from CuGeOs has been observed both, in the dimerized as 
well as in the uniform phase |@-^. In the low temperature uniform phase, for Tsp < T <^ Ji, 
the Raman spectrum displays a broad continuum centered at hu ~ 2.4Ji. Early on this 
continuum has been related to four-spinon excitations P-|^JT8[| . In the high temperature uni- 



form phase, and in addition to the four-spinon continuum, the spectra show a pronounced 
quasi-elastic line. In the dimerized phase the Raman intensity develops a gap at approxi- 
mately I.5-I.8A52- |2^ where A^^- is the singlet-triplet gap observed in INS |^. Moreover, 
four characteristic peaks appear in the spectrum the lowest one of which at 30cm~^ 'coin- 
cides' with the Raman gap. At present the interpretations of these peaks are controversial. 
Tentatively the 30cm~^ line has been attributed to a continuum of two-magnon bound states 
||19| , |20[| . A possible dimensional crossover effect has been invoked for the 225cm~^ line P,pO 



In this respect the influence of inter-chain coupling in CuGeOs is an open problem ^1 . 

Here I will focus on Raman scattering from the uniform phase of the J1-J2-S model, 
i.e. at 6 = 0. First I will describe a spinless-fermion mean-field theory to treat the J1-J2 
Hamiltonian. Next the Raman intensity is expressed in terms of a four-fermion correlation 
function. Finally results for the Raman spectra are compared with experimental findings 
and are contrasted against other theoretical approaches. 



II. MEAN FIELD THEORY OF THE Ji J2 MODEL 

The spinless-fermion mean-field (MF) theory for the J1-J2 model is based on the 
Jordan-Wigner (JW) representation [^] of the spin algebra, i.e. S[ = (qq — ^) and 
S^ = (— )' nj</(l ^ ^j^j)^ii where q are fermion operators. Inserting this into (|I]) at 5 = 
one obtains 

+a(— qq_|_2q+iq+i — ^+1^+1^+2^; — CiCi + qqq_,_2q+2 + t)J , (2) 

where Ji is set to unity in section ^ and |T|. / runs over the lattice sites. The terms 
proportional to a are absent in the JW representation of the isotropic n.n. Heisenberg 
model. In contrast to the n.n. model, both, the transverse as well as the longitudinal n.n.n. 
exchange-interaction lead to four-fermion vertices. Longer range spin-exchange leads to 
even higher order couplings, as is obvious from the JW representation which implies 21- 
fermion vertices for /-th.-nearest-neighbor spin-exchange. To proceed I treat Hamiltonian 
(0) in MF-approximation. Allowing for all contractions of type (qc^) one gets 

Hmf = E{-[1 + (^ + B){1 - 2a)] cos(A;) 4c, 

k 

+ ^(^-5)(l + 2a) sin(A;)(cK+^-cI+,c,)} + const. , (3) 

where k is the momentum, A = A* = {c2iC2i^i), and B = B* = (cg^+j^Cgi). In principle 
contractions of type D = D* = {(^21+2^21) "^^ occur, however, their selfconsistent value can 
be shown to vanish identically. Note that (^ allows for both, a uniform, if A — 5 = 0, 
and a gaped phase, ii A — B ^ 0. Here no attempt will be made to describe the spin- 
dimer state using A ^ B and the relevance of the gaped solution of the MF-theory will 
be discussed elsewhere. In the uniform case (^ resembles a single-band spinon gas with a 
hopping amphtude t{T, a) to be determined selfconsistently 

Hmf = J2 ^kclck 



ef, = -t{T, a) cos{k) = -[1 + 2A{1 - 2a)] cos(fc) (4) 

^ = ^ E cos(A;)/(e,) , 

0<fc<7r 

where /(e) = [exp(e/T) + 1]^^ is the Fermi function. At zero temperature the spinon 
dispersion simphfies to 

e;.(T = 0) = - [1 + 2(1 - 2a) /n] cos{k) . (5) 



For vanishing frustration this is identical to Bulaevskii's result —(1 + 2/tt) cos{k) [Q. The 
latter is known to compare reasonably well with the exact spinon dispersion e^ = —7i/2 cos{k) 
by des Cloizeaux and Pearson [^. For finite a the MF-theory results in a frustration 
induced softening of the spinon stiffness which, for T = 0, can be expressed as 

V (a) 4 

^V^ = 1 a ^ 1 - 0.778a , (6) 

Vs{0) 2 + 7r ' ^ ^ 

where Vs{a) = dek/dk\k=n/2 is the spinon velocity. Eqn. (|^) is qualitatively consistent 
with a numerical study where Vs{a)/vs{0) ~ 1 — 1.12a has been found fl^. In fig. |l| 



t(T, a) is depicted as a function of temperature for various values of frustration. As is 
obvious the MF spinon-stiffness is decreased, both, as a function of increasing frustration 
and temperature. At zero temperature the MF ground state energy is given by EqIo) = 
N{-l/7i - I/tt^ + 2a/7r2) which, for a = 0, leads to Eo{0) ^ -0.420 [El. This agrees 



reasonably well with the Bethe-Ansatz result Eq = 1/4 — ln(2) ^ —0.443. Moreover, the 
linear frustration dependence of the MF ground state energy, (Eo^a) — Eq{0))/N ^ 0.203a 



is close to that found in finite-chain diagonalization p5[, where {Eo{a) — Eq(0)) / N ^ 0.177a 
for a < 0.3. 

III. RAMAN SCATTERING 

The MF-theory is a convenient tool to study Raman scattering from the J1-J2 model. 
The Raman vertex is given by Loudon-Fleury's photon-induced super-exchange operator 

m 



-R — 2^ Tlmi^in ■ n;m)(Eout " n.im)Si ■ S^ • (7) 

Im 

Here Tim sets the coupling strength, Ej„ (Eout) refers to the field of the in(out)going light, 
and nim labels a unit vector connecting the sites / and m. By symmetry, in a strictly ID 
situation, (|^) leads to scattering only for parallel polarization along the c-directed chains. 
Since the real-space decay of Tim is comparable to that of the exchange integrals in (|I]) it is 
sufficient to consider a Raman operator with at most n.n.n. spin-exchange 

i? = C^(Sz-Sz+i+/?S;-Sz+2) . (8) 

I 

The scattering intensity I{uj) is obtained via the fiuctuation dissipation theorem I{uj) = 

x"{uj)/(l — e~^^'^) from the dynamical susceptibility of the Raman operator 



X"(^) 



r+oo . 

ImMuj + iri)] =Re dt e'^^^'"^' {[Rit) , R]) , (9) 

JO 



where (...) denotes the thermal average. 

Frustration of unequal magnitude regarding the J1-J2 model and the Raman operator is 
mandatory for non-vanishing inelastic scattering. li a = (3 the Hamiltonian for 6 = and 
the Raman operator commute which leads to elastic scattering only. The inelastic intensity 
which results from (H) is identical to that of a Raman operator R' = R — 7/J. Setting 
7 = C(3/a or 7 = C leads to scattering by a renormalized n.n. or n.n.n. Raman operator 
only 

Ri = C{l-^)Y.Si-Si+r or R^ = C{(3 - a)Y^Si -Si^^ • (10) 

« T T 

At present, an exact treatment of (j^) is not feasible and the Raman intensity of any ap- 
proximate evaluation will depend on 7. In this respect R2 is the proper choice for a, /? -C 1 
since it guarantees that /(cij) oc (/3 — a)^. This is not obvious for Ri or other values of 7. 
In those cases, even for a,/? ^ 0, R' may contain a component of order unity, proportional 
to the Hamiltonian, which has to be projected out. To avoid this complication R2 will be 
considered hereafter. 
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Using the JW-ferinions the Raman operator R2 can be expressed in terms of a four- 
fermion operator 

R2 = C{P -a) J2 Kk, k', q)clck+gcl,Ck,_g + Aipt (11) 

k,k',q 

h{k, k', q) = cos(2g) — cos(2/c + q) — cos(2A;' — q) , 

where Aiph labels all one-particle-hole excitations of the JW-transform of /?2- Since these 
occur at zero total momentum they do not contribute to inelastic Raman scattering within 
MF-theory [^. Using (pTj), the Raman intensity (P) is written in terms of the four-fermion 
propagator 

X{r) = C\a - (3)' ^ [h{k, k\ q)h{p, p', r) 

k,k',q,p,p',r 

(T^(4(r)Cfc+g(r)4,(r)c;.,_g(r)c],,_^Cp,4+^Cp))] . (12) 

Here r is the imaginary time and x"{^) results from the usual analytic continuation x"{^) = 
—Im[x{iuJn -^ uj + iO^)] where Un = 2n7TT is a Bose Matsubara-frequency. On the level of 



MF-theory the four-fermion propagator ([T^ ) is evaluated neglecting all vertex corrections 



and using the MF one-particle Green's functions corresponding to (§). After a number of 
standard manipulations I obtain 

1 r ^ r ^, ^\ g\k,k' + q,q) 






[2tsin(g/2)]2 - {ek+g - ek - uj) 



[/(cfc) - /(e/c+g)] [/(cfc') - /(cfc'+g)] [«(efc'+g " Cfc' + uj) - n{ek'+q - ek')] > , (13) 
where g{k, k', q) = \[h{k, k', q) + h{k, k', k' — k — q)] and the discrete sum on k' runs over 



2 
the set of solutions of 



Cfc+g — Cfc — ^ 



-'^'+«/^'-^Sh;(^ ' *'^i-.-i ■ (") 



This concludes the MF-theory of Raman scattering. 



IV. DISCUSSION 

In fig. 1 compare the MF-intensity obtained from a numerical integration of (0) (solid 
and dashed lines) at T = 20K and 12.8K, for Ji = 160K and two values of frustration 
a = 0.24 and a = 0.35 with an experimental Raman spectrum observed at T = 20K in the 
uniform phase of CuGeOs |0,|TB[ (dashed-dotted line with markers). Phonon lines at 184cm~^ 
and 330cm~^ have been removed. This Raman data is consistent with other published results 
0J^. The absolute magnitude of the observed Raman intensity as well as the coupling 
strength C of the scattering operator (^ are unknown quantities. Therefore the two sets 
of MF-intensities for distinct a and the experimental spectrum have been normalized on 
a scale of arbitrary units. The theoretical parameters in fig. § represent no attempt at a 
'best fit' to the data, rather they have been chosen among those obtained independently 



from studies of the magnetic susceptibility of CuGeOs pOHT^] • Figure g displays reasonable 
agreement between theory and experiment, although details of the experimental spectrum, 
i.e. the structure at 375 wave numbers and the initial curvature are not reproduced by the 
MF-theory. Evidently a value of a = 0.24 |T^ is slightly more favored by this comparison 



at Ji = 160K than a = 0.35 pT| , p!2| . Similar findings have been made in numerical studies 
pO| . However, since the maximum of the MF-intensity roughly scales with the maximum 
of the spinon-dispersion oc Ji[l + 2(1 — 2a) /it], a better agreement with a = 0.35 instead 
of a = 0.24 can also be reached simply by increasing Ji by ~10%. 

The comparison in fig. |^ strongly corroborates a study |18| of Raman scattering from 



the uniform phase of CuGeOa which is based the 'solitonic' mean-field description of the 
XXZ Heisenberg chain by Gomez-Santos |2^. This approach employs a domain-wall rep- 
resentation of the spin chain which leads to formal developments very different from the 
straightforward application of the Jordan- Wigner type of MF-theory presented here. 

In fig. 1^ the temperature dependence of the MF-spectrum is shown, both, in terms of the 
intensity I{uj) - which is observed in experiment - and the Raman-operator susceptibility 
x"{uj). The latter exhibits a left-shift of its maximum upon increase of the temperature. 

8 



This leads to enhanced low-frequency spectral weight which, by virtue of the Bose-prefactor, 
turns into a quasi-elastic line in I{uj) for T>Ji. A corresponding tendency has been de- 



tected in finite-temperature Lanczos studies |32|. In MF-theory the high-temperature en- 



hancement of the low-frequency spectral weight in x"{^) is due to the reduction of the 
spinon-stiffness as a function of increasing temperature, see fig. |I], and due to the 'smear- 
ing' induced by the Fermi- and Bose-functions contained in (0). It is tempting to relate 
the behavior depicted in fig. |^ to the high-temperature quasi-elastic line observed in the 
uniform phase of CuGeOs pj-Q. However, in CuGeOs the four-spinon continuum is only 
weakly shifted by temperature and remains more clearly separated from the quasi-elastic 
line for Tsp <T <Ji. This discrepancy may be due to an overestimation of the temperature 
dependence of the spinon-stiffness in MF-theory |29|j3^ . 

Finally, I emphasize that the MF-theory gives only a limited description of spinon in- 
teraction effects. In the case of the two-spinon propagator this is known to result in an 



incorrect description of the spectral weight distribution ||31| . This caveat of MF-theory has 



stimulated a study of approximate vertex corrections to the /our-spinon-propagator [52 



However, the resulting Raman spectra show no agreement with experiment. At, present 
the impact of spinon interaction effects beyond MF-theory on the Raman spectra remain 
unclear. 

In conclusion I have described a finite temperature MF-theory for frustrated spin- 
chains and consequently detailed its application to magnetic Raman scattering. At low- 
temperatures, in the gapless phase, I find a scattering continuum which is due to frustration 
induced four-spinon excitations and is compatible with observed Raman spectra of CuGeOs. 
In the high-temperature regime the MF Raman-intensity is dominated by a quasi-elastic 
line which results from the temperature dependence of the spinon spectrum. 
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FIGURES 
FIG. 1. Temperature dependence of the mean-field hopping amplitude for various values of 

frustration. 

FIG. 2. Mean-field Raman intensity (solid and dashed lines) at two temperatures and for 
two values of frustration as compared to the experimental Raman spectrum of CuGeOs [§,18| 
(dashed-dotted line & markers). 

FIG. 3. Mean-field Raman intensity and susceptibility for various temperatures. 
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